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ABSTRACT 


A   computer    simulation    study    of    channeling    in   a    diamond 
lattice.       The    simulation   was    done   for    a   xenon    ion    striking 
the    (100),     (110)    or    (111)    surface    of   a    silicon    target. 
Potential    functions    for    the   Si-Si    lattice   bond   and    the   Xe-Si 
interaction    are  postulated.      The    electronic    stopping    cross 
section    for    the   (110)    channel    of    silicon    is    estimated. 

This    work    is    a    continuation    in    the    development    of  a 
computer    model    formulated   at    the   USNPGS   which    takes    into 
consideration    the    displacement    of    the   atoms    in    the    target 
lattice  as    well    as    inelastic    energy    losses    by    the  primary 
ion.       The    lattice    was    not    thermalized   and    only    the   repulsive 
portion    of    the    lattice-lattice  potential   was    utilized. 
Computed   ranges    are    in    good   agreement    with   experimental    date 
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I .   INTRODUCTION 

Experimental  and  theoretical  studies  over  the  past  decade 
have  established  that  the  slowing  down  of  ions  in  crystals  is 
strongly  dependent  on  the  cr ys tallogr aphic  orientation  of  the 
lattice   with  respect  to  the  direction  of  the  incident  beam. 
[.1,2,3,4,5]   This  effect  is  known  as  channeling.   A  particle 
moving  along  one  of  the  major  axis  of  a  crystal  may  be  "steered" 
(i.e.  its  momentum  redirected  toward  the  center  of  the  channel) 
by  successive  gentle  collisions  with  the  "strings"  of  atoms  that 

comprise  the  "walls"  of  the  channel.   Davies  et  al  [ l]  reported 

125 
that  this  channeling  effect  increases  the  mean  range  of  Xe 

by  as  much  as  a  factor  of  10  over  the  values  observed  in  a 

non-crystalline  solid  of  similar  atomic  number. 

Earlier  computer  studies  performed  at  the  NPGS  investigated 
channeling  in  the  body  centered  [6]  and  face  centered  cubic 
lattice.  [7]   This  thesis  is  an  extension  of  this  previous  work 
into  the  diamond  lattice. 

When  a  substrate  is  bombarded  by  a  beam  of  incident  ener- 
getic particles  it  will  lose  some  of  its  own  ions  by  sputtering 
and  retain  some  of  the  incident  ions.   The  incident  ions  retained 
are  said  to  be  implanted,  and  the  technique  of  using  an  energetic 
ion  beam  to  introduce  ions  into  a  substrate  is  called  ion  implan- 
tation. [5]   This  doping  technique  in  silicon  should  be  very 
effective  due  to  the  extremely  open  nature  of  its  lattice  and 
its  large  channel  sizes. 
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The  difference  in  the  depth  of  penetration  of  energetic 
particles  along  the  principal  channels  in  silicon  may  be  attri- 
buted to  the  various  channel  sizes.   Figures  1-3  show  the  loca- 
tions of  the  various  channels  in  the  (100),  (110)  and  (111) 
surfaces  of  silicon  while  Fig.  4  is  a  "ball  and  stick"  model 
showing  the  relative  shape  and  sizes  of  the  three  principal 
channels  in  silicon.   According  to  Davies  [  l]  and  Manchester  [8] 
the  preferential  order  of  channeling  would  occur  along  the  ( 110) , 
(ill)  and  ( 100)  channels  in  decreasing  order  with  the  (111) 
channel  only  marginally  better  than  the  (100)  channel.   This  fact 
is  graphically  indicated  in  Fig.  5- 

Much  work  over  the  past  several  years  has  been  devoted  to 
studying  this  channeling  effect  in  silicon.   Manchester  [8]  has 
studied  the  feasibility  of  forming  metallurgical  junctions  in 
silicon  by  making  use  of  the  channeling  phenomenon.   Eisen  L  9~) 
has  studied  the  channeling  of  medium  mass  ions  through  silicon 
with  special  attention  directed  toward  the  electronic  stopping 
cross  section  of  well  channeled  ions.   Davies  and  his  co-workers 
have  written  a  series  of  papers  [2,3,4,10]  discussing  the  pene- 
tration of  KeV  and  MeV  projectiles  in  silicon  and  the  disordering 
of  the  lattice  caused  by  this  ion  bombardment.   Glotin  J.  11~]  , 
Nelson  and  Mazey  [ 12] ,  and  Dearnaley  et  al  [ 13]  have  investigated 
the  implantation  of  phosphorous  in  silicon.   Eriksson  et  al  f 14] 
have  studied  the  implantation  and  annealing  behavior  of  Group  III 
and  V  ions  in  silicon.   Gibbons  f5]  sums  up  the  current  theoreti- 
cal and  experimental  work  pertinent  to  the  problem  of  predicting 
impurity  distribution  profiles  in  implanted  material. 
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II .   STUDY  OBJECTIVES 

This  report  is  an  extension  of  the  NPGS  channeling  model  from 
the  bcc  and  fee  orientations  into  the  diamond  lattice.   This 
thesis  attempts  to  establish,  with  xenon  as  the  bullet  and  sili- 
con as  the  target,  the  following  characteristics  for  the  three 
major  channeling  orientations: 

a)  The  potential  between  the  bullet  and  the  lattice, 

b)  A  value  for  the  electronic  stopping  power,  and 

c)  The  ion  implantation  profile  at  5  KeV . 
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III.   SIMULATION  MODEL 

The  basic  model  has  been  developed  by  the  NPGS  group  over 
the  past  six  years.   A  brief  description  of  this  model  and  the 
alterations  made  for  this  study  follow  but  for  detailed  pre- 
sentations the  reader  is  referred  to  Refs.7,15  and  16 . 

The  simulation  model  consists  of  a  single  primary  and  a 
silicon  (diamond)  lattice  target  containing  between  60  and  150 
atoms.   All  runs  were  made  with  the  primary  (bullet)  striking 
the  target  normally,  although  the  program  is  sufficiently 
flexible  that  the  bullet  could  be  fired  into  the  target  at  any 
angle. 

The  silicon  target  lattice  as  used  in  this  simulation  is 
composed  of  a  central  core  and  a  surrounding  shell.   The  program 
is  structured  so  that  the  bullet  is  never  allowed  to  leave  this 
central  core.   In  addition,  in  the  (110)  and  (111)  orientations 
the  shell  atoms  were  not  permitted  to  move.   This  allows  the 
program  to  run  slightly  faster  and  also  tends  to  hold  the  core 
atoms  in  place  for  a  slightly  longer  time.   The  shell  atoms  were 
not  allowed  to  displace  in  the  (100)  channeling  direction.   Even 
with  the  shell  atoms  displacing,  given  a  bullet  energy  and  an 
equivalent  channeling  location  in  the  (111)  and  (100)  orien- 
tations, the  (100)  channeling  program  ran  faster.   For  a  detailed 
discussion  of  the  lattice  generator  to  include  the  dividing  of 
the  target  crystal  into  active  core  and  passive  shell  see 
Appendix  A. 
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In   previous    computer    simulations    the    interatomic   potential 
function   and   force   function   were   of    the   exponential    (Born-Mayer) 
type 

F   =    exp(A    +    Bx) 

where  A  and  B  are  empirically  determined  constants  and  x  is  the 
atomic  separation.   The  functions  are  not  always  of  this  type 
in  this  simulation,  and  force  and  potential  tables  were  con- 
structed.  For  details  of  this  scheme  see  Appendix  B. 

Two  major  assumptions  are  made  in  this  simulation  model, 

(1)  The  bullet  moves  in  a  perfect  lattice  undisturbed  by  the 
thermal  displacements  of  the  atoms . 

(2)  Only  repulsive  forces  between  the  lattice  atoms  have 
been  included. 

While  both  of  these  assumptions  are  obviously  erroneous, 
previous  studies  have  shown  that  this  simulation  model  is  still 
a  good  first  order  approximation  for  channeling. 
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IV.   PROCEDURE 

This  simulation  study  was  conducted  in  an  IBM  360/67  computer 
using  FORTRAN  IV  language. 

The  computer  sets  up  the  diamond  lattice  in  the  desired 
orientation  (see  Figs.  1-4)  and  starts  a  xenon  ion  into  it.   A 
brief  study  of  the  symmetry  of  the  target  lattice  shows  that  the 
indicated  impact  areas  cover  all  the  possible  points  for  the 
three  orientations  of  silicon.   Within  these  impact  areas, 
several  points,  depending  on  the  orientation  of  the  lattice,  were 
selected  as  representative  of  the  channel.   In  the  determination 
of  how  many  points  were  to  be  chosen  in  each  impact  area,  con- 
sideration was  given  to  obvious  lattice  symmetry  and  to  computer 
running  time. 

Various  combinations  of  Bullet-Target  Potentials  and  CELS, 
a  constant  used  to  determine  inelastic  energy  losses,  were  tried 
at  different  bullet  energies  and  the  maximum  range  of  the  bullet 
for  each  orientation  was  matched  against  experimental  range  data 
as  reported  by  Davies  et  al  [ l] . 

Channel  profiles  for  5  KeV  xenon  atoms  were  determined  for 
the  (100),  (110),  and  (111)  channels  of  silicon.   Runs  were  made 
to  determine  dE/dx  for  each  major  axis  in  silicon  and  maximum 
ranges  for  xenon  at  higher  energies  were  determined  knowing 

-/  <f  > 


1 
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V.       POTENTIAL    FUNCTIONS   AND    INELASTIC   LOSS    CONSTANT 

The   purpose    of    this    section       is    to   discuss    first,    the 
determination    of    the  potentials    used   to   approximate    the   Si-Si 
lattice    bond  and    the   xenon-silicon    bullet-lattice    interaction, 
and   second,    the    determination    of    an    inelastic    loss    constant. 
Presently,    there    exists    in    the    literature   calculated  potentials 
for    the   Si-Si    system  [  17]    but    none   for    the   Xe-Si    system. 
Eisen  I  9]    has   published  experimental   electronic   stopping   cross 
sections    for    medium   mass    ions    in    silicon    but    no    data    exists    for 
xenon    in    silicon. 

A.       Si-Si    BOND 

The    starting   point    for    selecting   a   potential    for    the   Si-Si 

bond   are   the    Interaction   Potential    Functions    (IPF)    as    calculated 

by   Harrison   [ 17] •      These  potentials    are   only    the   repulsive  portion 

of    the  potential    function    for    pairs    of  particles    and   neglect   any 

lattice    effects.       Harrison   [  17]    also    indicates    that    they    are   most 

likely    too    "hard".       Figure   6    shows    the   radial    electron    density 

for    Si.       Since    this    radial    electron    density    was    determined   using 

a    spherical    approximation,     it    is    not    correct    for    r   ^    1 . 5^ 

H 

(a       =    Bohr    Radius    =   O.529X)    [  17]  .       Figure   7    shows    several    dif- 
H 

ferent  potential    functions    for    the   Si-Si    bond   with    the  Gibson    II 
for    comparison.       The    Si-Si    and   Si-Si      IPF's    show   approximately 
the    correct    behavior    and  zero   close    to    the   nearest    neighbor    dis- 
tance   for    silicon.      The   Si       -Si    "    also    shows    approximately    the 
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correct    behavior    but    fails    to   zero   at    the  nearest   neighbor 

+1         +1 
distance.       Curve    3    of   Fig.    7    is    the    Si      -Si         IPF   forced   to 

zero   at    the   nearest    neighbor    distance. 

+  + 

Any    one    of    these    three  potentials    (Si-Si,    Si-Si      or    Si-Si 

(zeroed))  will  serve  as  a  good  approximation  for  the  Si-Si  bond. 
It  is  obvious  that  the  Gibson  II  will  not  serve  as  a  good  repre- 
sentation   of    this    bond. 

The   covalent    bonding    between    Silicon    atoms    has    the   effect 
of   concentrating    about    0.5    of   an    electron    of    the    outer    shell 
along    every    chemical    bond  ll83.      Therefore,    it   seems    reasonable 
to   conclude    that    each    silicon    atom   along    the   Si-Si    bond   sees    the 
other    in    a    slightly    ionized   state   but    not    completely    ionized 

to    the    +1    state    due    to    the   shielding    effects    of   this    shared 

.     .  .  .+1       .+1 

electron.      This    would   seem    to    eliminate    the   Si-Si    and    Si      -Si 

potentials    for    the   Si-Si   bond   and   lead    to    the   selection    of    the 

Si-Si      potential    as    representing    the   bond.       It   was    noted   earlier, 

however,    that    the    IPF's    are  probably    too   hard.       Therefore,    the 

Si      -Si         (zeroed)    potential   was    chosen   as   representative   of    the 

+ 
bond  potential    as    it    is    "softer"    than    the   Si-Si      potential    and 

is    forced    to   zero   at    the   nearest    neighbor    distance. 

.  +         + 
Admittedly    the   reasons    for    choosing    the   Si    -Si      potential 

were  rather    arbitrary    but    the  reader    is    asked   to   keep    three 

points    in    mind.       First,    all    three   potential  are   approximately 

equal    and    trial    runs    indicated   no    significant    effect    on    the   range 

when    these   potentials    were    interchanged.       Second,    the    neutral 

silicon    atoms    are    too    large   to   fit    into    the   Si    lattice.       Third, 
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the    two   overwhelming    factors    that    control    channeling   ranges    are 
the    bullet-lattice  potential    and    the    inelastic    loss    approximation. 
Any    reasonably   potential    selected   for    the   Si-Si    bond   should   have 
no    discernable   effect    on    these   ranges. 

B.   INELASTIC  LOSS  CONSTANT 

One  of  the  principal  mechanisms  of  energy  loss  for  channeled 
ions  in  a  solid  is  the  interaction  of  the  ion  with  electrons 
(both  bound  and  free)  of  the  solid.   This  computer  study  accounts 
for  this  loss  through  the  use  of  a  "frictional  force  multiplier" 
which  acts  on  the  ions  velocity  in  a  manner  to  produce  a  "drag" 
force  on  the  incident  ion.   This  "friction  force  multiplier"  is 
directly  proportional  to  the  electronic  stopping  cross  section. 

Eisen  L  9]  reports  electronic  stopping  cross  sections  (S  ) 

for  medium  mass  ions  (5  <  Z  ^  19)  in  the  < 100) ,  (110),  and  <111> 

channels  of  silicon.   In  unpublished  data  [  19]  he  extends  S   for 

e 

the  (ll())     channel    by   allowing    the   atomic    number    of    the    incident 
ion    to    increase   up    to   Z   =    3°    (Krypton).       Eisen's    data    reveals    the 
following:    first,    S      was    found   to    exhibit    a    strong    oscillatory 
dependence    on    Z,    the   atomic    number    of  the  incident    ion,    second,    S 
is    lowest    in    the   (110)    channel    and   third,    S      is    approximately 
equal    in    the   (ill)    and   (100)    channels. 

No   data   presently    exists    in    the    literature    for    the   electronic 
stopping    cross    section    of   xenon    in    silicon.      However,    based    on 
previous    computer    studies    of    channeling    in    copper    and    tungsten 
undertaken   at    NPGS    the   frictional    force   multiplier    was    estimated 
to    be    1.0   x    10       J    n-sec/m    in    the  ( 100)    channel.       Using    Eisen's 
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data  [19]  and  noting  the  ratio  between  S   in  the  (100)  and  (110) 

channels,  the  frictional  force  multiplier  in  the  (110)  channel 

is  approximately  -8  x  10     n-sec/m.   These  values  were  used  in 

this  computer  study. 

Ericksson  et  al  ["20]  give  estimates  of  electronic  stopping 

cross  section  for  various  ions  along  the  (100)  and  (110)  channels 

of  tungsten.   Their  data  also  indicates  that  S   shows  a  strong 

e 

oscillatory    dependence    on    Z,    the   atomic    number    of    the    incident 
ion    and   this    oscillatory    dependence    is    very    similar    to    the 
dependence   Eisen   [9,19]    found   for    ions    channeled    in    silicon.       If 
the   assumption    is    made    that    the   ratio   between    stopping    cross 
sections    for    xenon   and   krypton    (which    is    chemically    similar    to 
xenon)    channeled    in    tungsten    is    approximately    the    same   as    the 
ratio   between    these    ions    channeled   in    silicon,    one    obtains    a 
value    of 

8.5   x    lO"1^    eV"Cm         (where   v    =    1.0   x    10      cm/sec) 
atom         v  ' 

as    an    estimate    of   S      in    the   (110)    channel    in    silicon.      This 

e 

-14 
converts    to  a    value    of   7   x    10  n-sec/m   as    a    frictional    force 


mu 


ltiplier    for    the    (110)    channel    and    is    in   good   agreement   with 


the   value   used    in    this    computer    study. 

se  (v)  -TM.6   x    10~19(j/eV)         .  -  , 

CELS    =  ' — ' ; — Sr ~      where  T\    =    number    density. 

v(cm/sec)    x   10       (m   /cm    ) 

The   author    does    not   make   any    claim    that    the   values    of    the 
"frictional    force"    used    in    this    study   are    "correct".       It    is    felt, 
however,     that    based   on    a    study    of   Eisen's    data    [9519]    and    the 


extrapolation    described   above,    the    values    chosen    for    S      will 

e 

certainly    be   within    20%    of    the   experimentally    determined    values 
and   more    than    likely   will    be   much    closer    than    that. 

C.       Xe-Si    POTENTIAL 

As    in    the   case    of    the    Si-Si   bond,    no  potential    function 
exists    for    the   xenon-silicon    interaction.      Therefore,    the 
Interaction   Potential    Functions   will    again    serve   as    a    starting 
point    for    selecting    the   Xe-Si   potential.      Figure    8    shows    the 

radial    electron    density    for    xenon    in    the   nonionized  and   +6 

.       .         .  .  .  +6        +4 

ionization    state,    While   Fig.    9    shows    the   IPF   for    Xe       -Si        and 

+6       .+1 
Xe      -Si 

Far    from    the   Si-Si    bond,    the   electron    density    for    silicon 

+4 
should   approach    the   Si         distribution    (Fig.    6).       Special    note 

should   be    taken    of    the   (111)    channel    of   silicon   where    the   atoms 

bounding    the    channel    have    one    of    their    bonds    parallel    to    the 

centerline    of    the   channel,       Determing    the    ionization    state    of 

xenon    is    also    difficult.      Harrison    [ 17]    has    found   for    Xe-W 

system   a    good  potential    fit    can    be    found    if  xenon    is    assumed   to 

be    in    the    +6    or    +8    ionization    state,    even    though   this    high    ionic 

+      +6 
state    is    somewhat    surprising.      The   Xe    -W        potential    was    found 

to   be  unr ealistically    strong.       As    a    result,    of    the  above    it   was 

+6+4  . 

concluded    that    the   Xe       -Si         IPF   was    a    good   starting   point    for 

the   bullet -lattice   potential. 

+6       .+4  .  .  .       . 

The   Xe      -Si         IPF    (Fig.    9)    diverges    to   unr ealistically 

strong    values    at    large    separations.      This    divergence    is    charac- 
teristic   of    the    theoretical    model    used   to   calculate   the   IPF 
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and    is    not    characteristic    of    the  physical    system   [ 17] •      A   Born- 
Mayer    approximation    for    the   Xe-Si    interaction   would   fail    at    small 
separations    because    it   would   approach   a    finite   value   at    zero 
separation,    while    the   actual    Xe-Si    function    should   approach   a 
pure   coulomb  potential    at   small    separations.      Our    final    choice 
of   the  potential    to    represent    the   Xe-Si    system    incorporated   both 

of    the   above   approximations.       For    separations    less    than    1.2a 

H 

the   IPF   was    used,    while   a    Born-Mayer    approximation   was    constructed 

for    separations    greater    than    1.2a    . 

H 

Curve    BM1    (Fig.    9)    when    used    in    conjunction   with    the   pre- 
viously   discussed  values    of   the    "frictional    force"  parameter 
gives    excellent   results    for    xenon    ions    injected    in    the   (  110)    and 
(100)    channels,    however,     it    is    a   poor    approximation    for    xenon 
in    the   (ill)    channel.       After    several    trial    runs    it    became    obvious 
that   all    three    channeling    directions    could   not    fit    on    one   Born- 
Mayer    approximation.       The   maximum    range    in    the    (ill)    channel 
could   be   forced    to   approach    the   experimental    range    if   the 
"frictional    force"    constant    is    decreased   by    two    orders    of   magni- 
tude.      This    is    unrealistic    since    it    conflicts    drastically  with 
Eisen's    [9]    experimental    data. 

A   solution    to    the  problem  was    to   make   a    separate   Born-Mayer 
approximation    for    the   (111)    channel;     this    is    curve   BM2    in    Fig.    9. 
Intuitively,    one    feels    that    this    answer    is    not    the   correct    one. 
All   previous    studies    have    found   that    one  potential    satisfies    all 
channeling    directions.       However,    two   points    should   be   kept    in 
mind.       First,    all   previous    computer    studies    dealt    with   channeling 
in    a    metallic    target   where    the  principal   method    of   bonding    between 
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lattice   atoms    is    metallic.      This    is    the   first    computer    study   where 
the    bonding    between    the    lattice   atoms    is    co-valent    and   any    con- 
clusions   drawn    for    a   metallic    target    need   not    necessarily    be    true 
in    silicon.       Secondly,    the    orientation    of    the   bond   between    silicon 
atoms    with    relation    to    the   channeling    directions    varies    among   all 
three  principal    channels    of    silicon.      As    previously   mentioned, 
the    extreme    is    reached    in    the   (ill)    channel    where    the    Si-Si    bond 
is   parallel    to   the    center    line   of    the   channel.      This    more   than 
likely    has    an   effect    on    the   silicon    charge    distribution    as    seen 
by    the   xenon    ion    and    could  possibly    cause    the   potential    for    the 
(111)    channel    to    differ    from    the    other    channels. 

Another    possibility    is    that    a    Born-Mayer    approximation    is    not 

+6+4  . 

valid    in    silicon.       A  picture   of    the    Xe      -Si         potential    might   be 

obtained   by    combining    the  portion    of    the   three    curves:     IPF,    BMl 

and   BM2    that    we    know   give   us    accurate    results.       Such   a    composite 

curve    is    shown    in    Fig.    10.       It  is  obvious    that    this    composite   curve 

is    unrealistic. 
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VI.       RESULTS 

A.       (100)    CHANNEL 

eV 
Figure    11    is    a   plot    of   dE/dx( in      /LU)    versus    ion    energy    for 

a   perfectly    channeled   xenon    ion    in    the   (100)     direction.       (Note 

1LU   =   2. 71$    =   0.0646  ^g/cm2) .       Between    30-100   KeV,    dE/dx    is 

approximately   proportional    to    the    energy    indicating    that    in    this 

region    the    elastic    and    inelastic    collision   processes    are   both 

exerting    a    significant    influence.      This    agrees    with   experimental 

data    reported   by    Da vies    [  l]  . 

Figure    12    is    a   plot    of    dE/dx    versus    E2.       Above   800  KeV 

h    ■       ■ 
dE/dx    is   proportional    to    E      indicating    that    as    expected   the    in- 
elastic process    is    the  primary    loss    mechanism.       Both   Figs.    11 
and    12    show    that    below   20   KeV    the    elastic    collision   process    is 
primarily    responsible    for    slowing    down    the    ion. 

Figures    13    and   14   are   range  profiles    and   range  profile 
sections    for    the   (100)    channel.      Figure    13    shows    only    one-half 
of    the    channel,    the   range   profile   for    the    other    half    is    very 
similar    but    not    identical    to    this. 

Table   1    (below)    is    a    comparison    between    experimental    ranges 

obtained   by    Davies   [  l]    and   computer    calculated   ranges    for    the 

(lOO)    channel.      The   column    labeled   dE/dx    is    the   range    value 

obtained   by    using    the    dE/dx    curve    in    Fig.    11    and   the  relationship 

E 

0 


=j       (dE/dx)"1    dE  • 
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The    computer    calculated   range   at    40  KeV ,    although   higher    than 
the  exper  imental    value    reported   by    Davies,    falls    within    his    +3% 
error    estimate. 

Range    ( 100) 


Energy 

Exper  imental 

dE/dx 

Computer 

5 

60 

69.7 

20 

491.7* 

485 

490.2 

40 

1162.7 

1230 

1193-8 

Table  I.   Energy  in  KeV,  range  in  LU 

*Exper imental  range  for  20  KeV  extrapolated 
from  Davies  [4]  experimental  data. 

In  a  bcc  or  fee  crystal  an  ion  injected  into  the  center  of 
the  channel  remains  there.   This  is  not  the  case  for  an  ion 
injected  into  the  center  of  the  (100)  channel  of  silicon.   The 
ion  "wobbles"  slightly  around  the  channel  center  oscillating 
between  (110)  planes.   (Fig.  15)   Note  the  similarity  of  the 
orbit  in  the  second  and  fourth  quadrant. 

Figure  6  is  a  Range  Profile  Section  for  a  5  KeV  silicon  ion 
injected  into  the  (100)  channel  of  silicon.   Although  no  experi- 
mental data  exists  to  compare  maximum  ranges  it  is  interesting  to 
note  that  the  range  profile  section  is  very  similar  to  the  range 
profile  of  xenon  in  silicon. 

B.   (110)  CHANNEL 

Figure    7    is    a   plot    of   dE/dx    versus    energy    at    various    values 
of   CELS    (analogous    to    electronic    stopping  power)    for    a    xenon    ion 
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injected  into  the  center  of  the  (llO>  channel.   It  is  apparent 
that  the  energy  loss  is  very  sensitive  to  CELS  especially  at 
higher  energies  indicating  that  electronic  stopping  is  the 
principle  mechanism  of  energy  loss  for  well  channeled  ions  in 

the  (110)  direction.   Figure  l8  is  a  plot  of  dE/dx  versus 

h  -lk 

(ion  energy)   for  well  channeled  ions  with  CELS  =  -  8  x  10 

n-sec/m,  the  value  selected  for  use  in  this  study.   dE/dx  is 

proportional  to  E2  at  incident  ion  energies  greater  than  2  KeV 

again  indicating  that  electronic  stopping  is  the  principle  loss 

mechanism . 

Figure  19  and  Fig.  20  are  a  range  profile  and  a  range  profile 
section  for  the  (110)  channel.   Both  show  only  half  the  channel 
although  the  other  half  of  the  channel  should  be  identical. 
While  it  is  apparent  that  the  (110)  channel  is  a  preferred 
channeling  direction  the  maximum  range  falls  off  rapidly  if  the 
ion  is  injected  slightly  off  center.   (For  example,  if  the  ion 
is  injected  approximately  0 . 2LU  from  the  center  of  the  channel 
the  maximum  range  decreased  by  one-half).   Of  particular  interest 
is  point  1  in  the  Range  Profile  Section  (Fig.  20).   An  ion  in- 
jected at  this  point  in  the  channel  has  a  range  of  100LU  which 
is  much  greater  than  expected.   The  xenon  ion  injected  here  is 
strongly  influenced  by  the  silicon  atom  located  at  position  C 
which  forces  it  into  the  center  of  the  channel  and  adds  signi- 
ficantly to  its  range. 

Table  II  is  again  a  comparison  between  ranges  obtained 
experimentally  by  Davies  [  l]  and  computer  calculated  ranges. 
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Range  110 
Energy      Exper  imental      dE/dx       Computer 

5  --  1090       1094 

20  2371  2400 

Z+0  3812  38OO 

Table  II.   Ranges  in  LU,  Energies  in  KeV. 

The  ranges  calculated  using  dE/dx  are  in  good  agreement  with 
experimental  ranges  obtained  by  Davies  [ l] .   No  attempt  was  made 
to  channel  a  20  KeV  xenon  ion  as  it  would  be  too  expensive  in 
computer  running  time  (approximately  90  minutes). 

Figure  21  shows,  the  ranges  obtained  for  selected  points 
in  the  (110)  channel  when  the  energy  of  the  incident  xenon  ion 
was  increased  to  20  KeV.   Two  points  should  be  noted.   First, 
the  channel  "expands"  at  higher  energy  as  expected.   Second,  a 
comparison  of  the  range  of  the  ion  incident  on  point  1  with  the 
range  obtained  when  the  incident  ion  energy  was  5  KeV  (Fig.  20) 
shows  that  the  large  range  for  the  5  KeV  ion  is  a  low  energy 
phenomenon . 

C.   <1H)  CHANNEL 

Figure  22  is  an  enlargement  of  the  (ill)  channel.   Initially, 
point  1  was  thought  to  be  the  center  of  the  channel  but  subse- 
quent probing  of  the  channel  proved  this  assumption  to  be  incor - 
rect.   The  points  numbered  2  and  3  ar©  actually  the  channel 
"centers".   If  atoms  A,  B,  D  and  B,  C,  D  are  thought  to  form 
triangles j  as  indicated,  then  points  2  and  3  TOay  be  considered 
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the  center  of  their  respective  triangle  equidistant  from  the  atoms 
that  form  the  triangle's  vertices.   Additionally,  although  the 
basic  shape  of  the  channel  remains  unchanged,  the  atoms  that 
bound  the  channel  do  not  necessarily  lie  in  the  planes  indicated 
in  Fig.  22.   Two  other  possible  configurations  are:  (a)  atoms  A 
and  C  in  the  front  plane,  B  in  the  third  plane  and  D  in  the  fifth 
plane  or  (b)  atoms  A  and  C  in  the  fifth  plane  back,  D  in  the  third 
plane  back  and  B  in  the  front  plane.   These  three  channel  con- 
figurations are  shown  in  Fig.  3.   They  affect  the  ion  only  as  it 
initially  enters  the  lattice,  by  determining  in  which  direction 
it  will  be  initially  steered  and  have  almost  no  effect  (^  0.5%) 
on  the  calculated  maximum  range. 

Figure  23  and  24  are  plots  of  dE/dx  vs  ion  energy,  E  and 
E2.   Up  to  100  KeV  both  elastic  and  inelastic  collisions  are 
both  important  loss  mechanisms.   Below  20  KeV,  elastic  losses 
appear  to  predominate. 

Figure  25  shows  dx/dE  curves  for  xenon  in  the  (111)  channel. 
The  area  under  these  curves  should  be  an  accurate  representation 
of  the  total  range  travelled  by  an  incident  xenon  ion.   The  solid 
curve  (I)  was  obtained  by  injecting  xenon  ions  into  the  center  of 
the  (111)  channel  at  various  energies,  allowing  them  to  penetrate 
the  lattice  for  a  short  distance  (^  5LU)  and  determining  the 
energy  loss  of  the  xenon  ion  per  lattice  unit.   The  dotted  curve 
(II)  is  actually  a  20  KeV  run  for  a  xenon  ion  initially  injected 
in  the  center  of  the  (111)  channel  and  followed  until  it  finally 
stops  in  the  crystal.   The  difference  between  these  two  curves 
can  be  accounted  for  by  noting  that  an  ion  injected  into  the 
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center  of  the  (111)  channel  deviates  considerably  from  the 
channel  center  (up  to  .09LU)  as  it  passes  thru  the  crystal  and 
that  this  deviation  increases  as  the  ion  energy  decreases. 
Therefore,  the  implicit  assumption  made  in  determining  curve  I, 
that  the  ion  remains  in  the  center  of  the  (ill)  channel  as  it 
passes  through  the  crystal  is  incorrect  and  will  introduce  large 
errors  (up  to  20%)  in  ranges  at  the  low  end  of  the  energy  spec- 
trum where  the  elastic  loss  mechanism  predominates.   Table  III 
is  a  summary  of  maximum  ranges  obtained  at  various  energies  for 
the  (111)  channel. 

Range  <lll) 

Energy      Experimental      Curve  I  Curve  II     Computer 

5           --           37  60  57 

20         511         417  500  5H 

40        1209        1050  1170 

Table  III.   Energy  in  KeV,  range  in  LU. 

Figure  26  is  a  representation  of  the  (111)  channel  and  shows 
ranges  obtained  at  different  initial  impact  points  in  the  channel 
Two  things  can  be  noted  about  the  range  distribution.   First, 
there  is  a  definite  peak  at  the  center  of  the  channel  and 
secondly,  the  total  range  is  extremely  dependent  on  exact  im- 
pact point  in  the  channel.   This  last  fact  makes  it  impossible 
to  draw  any  meaningful  range  profiles  or  ranges  profile  section. 
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VII.   CONCLUSIONS 

1.  Maximum  channeled  ranges  were  fit  to  experimental  data 
at  20  KeV. 

2.  For  the  (100)  and  (ill)  channels  both  elastic  and 
inelastic  losses  are  important  between  20  and  80  KeV. 

3-   Inelastic  losses  become  predominant  at  energies  greater 
than  400  KeV  for  the  < 100>  channel. 

4.  Inelastic  loss  seems  to  be  the  predominant  loss  mechanism 
for  well  channeled  ions  in  the  (110)  channel  even  at  low  energies 
(~2  KeV) . 

5.  The  range  in  the  (ill)  channel  is  extremely  sensitive  to 
impact  parameter  and  appears  to  possess  no  symmetry. 

,       2 

6.  S   was  determined  to  be  8.5  x  10_1^  eV"Cm   in  the  (110) 

e  atom  ' 

channel    of   silicon.       In    the   ( 100)    and   (ill)    channels 

2 

^  -i  ^    ^         ,A-l4   eV-cm       .  ,     _.         ,  _.8  ,         . 

S      =    10.6    x    10  — — - —   (v    =    1.0  x    10      cm/sec). 

e  atom      v  ' 

7.  All    three   channeling    directions    for    xenon    in    silicon    can 
not    be    determined  using    the   same   Born-Mayer    Potential    Function. 
A   minimum    of    two    Born-Mayer    Potentials    are   required. 

8.  At    low    energies    the   (100)    channel    is    the   preferential 
direction    for    channeling   when    compared    to    the  (ill)    channel. 
This    is    the   reverse    of    the    situation    for    incident    ions    energies 
greater    than    20   KeV. 
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9.       Although    the   potentials    and    the    "frictional    force"    con- 
stant  were   discussed   separately    it    can   not    be   overemphasized 
that    all    three    of    these   variables    roust    be    treated   as    an    inte- 
grated whole.      A   change    in    either    the   bullet-target   potential 
or    the    frictional    force   constant    has    a    profound   effect    on    the 
maximum    channeled   range    especially    in    the   ( 110)    orientations. 
The   final    selection    of    the   Born-Mayer    approximations    to    the 
Xe-Si    interactions,    the   frictional    force  parameter    and   Si-Si 
bond  potential    seems    to   represent    the  best    of   all    the  parameters 
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APPENDIX  A  (LATTICE  GENERATOR) 

Previous  computer  studies  undertaken  by  the  NPGS  Group  over 
the  past  several  years  have  dealt  exclusively  with  the  body 
centered  and  face  centered  cubic  lattices.   In  building  the 
diamond  lattice,  which  is  basically  two  interpenetrating  face 
centered  cubic  lattices,  the  primary  aim  was  to  construct  the 
lattice  in  such  a  way  that  only  mimimum  modifications  on  the 
basic  computer  model  for  channeling  would  be  necessary.   The 
diamond  lattice  generator  was  constructed  so  that  no  major  change 
in  the  basic  lattice  regenerator  was  necessary.   Therefore,  to 
understand  the  lattice  generator,  it  is  necessary  to  understand 
the  fundamental  mechanism  behind  the  regenerator  subroutine. 

Subroutine  REGEN  works  as  follows: 

(1)  It  receives  the  parameter  MSS  from  the  main  program 
which  indicates  the  direction  of  regeneration  (+X,  +Y,  +Z) . 

(2)  It  checks  every  atom  in  the  crystal  against  a  constant 
which  depends  on  the  direction  of  regeneration  and  either 

a.  discards  the  atom  since  the  bullet  has  passed  and  is  no 
longer  under  its  influence,  or 

b.  shifts  and  renumbers  the  atom  while  preserving  its 
velocity  components. 

(3)  It  builds  and  undisturbed  section  crystal  in  the 
direction  of  travel  of  the  bullet. 

(4)  It  shifts  the  bullet  to  its  proper  location  in  the 
rebuilt  crystal  preserving  its  velocity  components. 
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Schematically,    this    is    done   as    follows:     (The   crystal    regen 
erates    when    the   bullet    passes   plane    E) 
A)    Crystal    Prior    to   Regeneration 
Direction    of  Travel 
Bullet 


A 


BCD 


y 


F      G      H 


Disordered  Crystal, 
Passed   By    Bullet 


The    follwing   planes 
are    identical: 

A    -    E    -    I 

B    -    F 

C    -    G 

D    -    H 


B)    Crystal    During    Regeneration 


_/v 


A 


BCD 


F      G      H 


Plane    E   goes    to   A 


F    " 

it 

B 

G   " 

it 

C 

H    " 

it 

D 

I    M 

it 

E 

Old  planes 

A,B. 

,c, 

D, 

discarded 

Discar  ded 


C)    Crystal    After    Regeneration 


New 
Plane 


Was 
Plane 


A 


/ 


Bullet 

B      C      D      E 


F      G      H 


F      G      H 


Planes  F,G,H,I  are  newly 
constructed  as  part  of  a 
perfect    crystal    lattice. 


This    scheme    discards    the  portion    of    the   crystal    the   bullet 
has    passed   through   and   builds    an    undisturbed  portion    of    the 
crystal    in    front    of    the   bullet. 
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To  build  a  crystal  which  will  fit  into  the  regeneration 
scheme  the  crystal  generator  must: 

(1)  Construct  one  plane  at  a  time  in  the  Y  direction.   It 
places  the  first  atom  in  the  Y  =  0  plane  at  the  lowest  value  of 
X  and  Z  (0,0)  and  continues  placing  atoms  at  increasing  values 
of  X  for  fixed  Z,  then  increases  the  value  of  Z  and  repeats  the 
process  until  the  plane  is  completed.   The  next  Y  plane  is  con- 
structed in  a  similar  manner  and  additional  Y  planes  are 
constructed  until  the  micr ocr ys tallite  is  completed.   It  assigns 
a  number  to  each  atom.   A  completed  Y  =  0  plane  for  the  (100) 

orientation  has  atoms  positioned  and  numbered  as  follows 

Z  fl2      *13     »14 

•  10     »11 

7      #8     ♦a 

•5     »6 


■3 J*> 


x 


Planes  must  be  generated  in  this  manner . 

(2)  Build  the  following  symmetry  into  the  crystal.   The  left 
and  right  face  are  identical  to  each  other  and  identical  to  one 
plane  in  between  them.   There  is  similar  symmetry  between  the 
front  and  back  faces  as  well  as  the  bottom  and  top  faces. 

This  symmetry  in  the  lattice  allows  the  regenerator  to 
operate  as  follows: 

(1)  It  checks  the  location  of  each  atom  in  the  crystal  to 
see  if  it  is  to  be  regenerated  (shifted). 
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(2)  If    it    is    to   be   regenerated,    a    fixed   number    (IDX,    IDY, 

or    IDZ    depending    on    the    direction    of   regeneration)     is    subtracted 
from   the   atom   number.      This    partially    shifts    the   atom    into    its 
identical    plane. 

(3)  The   regenerator    then    subtracts    (or    adds)    DXT ,    DYT ,    or 

DZT    from    one    of    theatom's    co-ordinates    depending    on    the    direction 
of   regeneration.      This    completes    the   movement    of    the  atom    into 
its    identical    plane. 

(4)  As   previously   mentioned   the    regenerator    preserves    the 
velocity   components    and  potential    energies    of   all    atoms 

r  egener ated. 

(5)  The  bullet  is  regenerated  in  a  similar  (but  not  identical) 
manner . 

In  addition  to  being  tailored  to  fit  SUBROUTINE  REGEN  the 
lattice  generator,  through  the  use  of  a  variable  LCUT ,  breaks 
the  lattice  into  an  active  core  and  a  surrounding  shell.   This 
is  an  absolute  necessity  due  to  the  large  size  of  the  crystal 
which  was  dictated  by  the  symmetry  requirements  of  SUBROUTINE 
REGEN.   (Orientation  (100)  contains  96  atoms,  orientation  (110) 
contains  60  atoms  and  orientation  (111)  contains  146  atoms). 
LCUT  does  this  in  the  following  way: 

1)  Any  atom  in  the  active  portion  of  the  crystal  is  given 
an  LCUT  =  0.   This  allows  it  to  interact  with  any  other  atom 
in  the  crystal . 

2)  Any  atom  in  shell  is  given  an  LCUT  =  1  which  allows  it  to 
interact  only  with  core  atoms  and  not  with  any  shell  atoms. 
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The  lattice  generator  subroutine  also  contains  the  constants 
XLL,  XLF,  ZLL;  ZLF  which  are  used  to  determine  when  the  lattice 
is  to  be  regenerated.   They  are  chosen  to  insure  that  the  bullet 
always  remains  in  the  core  and  the  lattice  regenerates  when  the 

bullet  is  about  to  leave  the  core  and  enter  the  shell. 

125 
As  an  example,  a  20  KeV  Xe     atom  in  the  center  of  the  (100) 

channel  of  silicon  has  a  range  of  491.7  LU  and  the  program  runs 

for  35  minutes  and  39  sees.   With  the  LCUT  package  inserted  into 

the  program  the  calculated  range  is  490.2  LU  with  a  program 

running  time  of  13  minutes  and  12  seconds.   The  LCUT  package  is 

obviously  effective  for  it  conserves  computer  time  without 

sacrificing  any  accuracy  in  range  calculations. 

In  summary,  the  Lattice  Generator 

(1)  Contains  all  the  constants  necessary  to  regenerate  the 
lattice, 

(2)  It  constructed,  with  appropriate  tests,  to  insure  that 
it  fits  into  the  scheme  necessary  for  Subroutine  REGEN,  and 

(3)  Contains  the  LCUT  package  which  reduces  the  number  of 
active  atoms  in  the  crystal,  thus  conserving  running  time. 

An  alternative  to  building  the  lattice  generator  to  fit 
subroutine  REGEN  would  be  to  construct  a  smaller  active  crystal, 
eliminate  the  LCUT  package  and  rewrite  the  REGEN  subroutine  to 
fit  the  lattice  generator .   This  method  was  not  chosen  for  the 
following  reasons: 

(1)  It  was  easier  to  modify  the  existing  lattice  generator 
to  fit  the  old  Subroutine  REGEN  than  to  start  from  scratch  and 
build  a  new  Regeneration  Subroutine  and  Lattice  Generator  Sub- 
routine tailored  for  each  other,  and 
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(2)  Simple  modifications  in  the  generator  alone  would  not 
drastically  alter  the  existing  Channeling  Program  which  had 
proved  to  be  an  accurate  representation  of  the  channeling  pro- 
cess in  fee  and  bee  crystals. 
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APPENDIX  B  (CALCULATION  OF  FORCES  AND  POTENTIALS) 


.on 


A  major  change  in  procedure  was  introduced  in  the  calculate 
forces  and  potentials.   In  previous  models  the  interatomic  po- 
tential function  and  force  functions  were  of  the  exponential  type, 

F  =  exp  (A  +  Bx) 
where  A  and  B  were  empirically  determined  constants  and  x  is 
the  atomic  separation.   The  atomic  separation  was  determined, 
substituted  into  the  appropriate  function  and  the  force  and 
potential  determined. 

This  computer  simulation  determined  forces  and  potentials  as 

follows.   First,  a  set  of  potentials  (in  Rydbergs)  derived  from 

Interaction  Potential  Functions  are  read  into  the  program  for 

the  lattice-lattice  and  bullet -lattice  interactions.   These 

potentials    start   at    0.01a       (a      =    Bohr    radius    =   0.529A)    and  are 

H         H 

spaced   every    0.01a      up    to   5»0a       interatomic    separation.      Secondly, 

H  H 

a  control  card  associated  with  each  deck  of  potentials  is  read. 
This  control  card  allows  the  extrapolation  of  a  Born-Mayer  type 
potential  (if  desired)  from  any  point  on  the  Interaction  Potential 
Function  Curve.  Thirdly,  the  potential  at  the  nearest  neighbor 
distance  is  determined  and  then  the  potential  is  zeroed  at  that 
distance.  All  potentials  at  atomic  separations  smaller  than  the 
nearest  neighbor  distance  are  adjusted  accordingly.  All  poten- 
tials at  separations  larger  than  the  nearest  neighbor  distances 
were    set    equal    to   zero.      This    does    not    leave   any   potential    free 
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regions  in  the  lattice.   Fourthly,  potentials  are  converted  to 

electron  volts.   Finally,  the  relation  between  the  force  and 

potential  energy  can  be  expressed  simply  as 

F(X)  =  -  VU(X)  . 

Assuming  that  the  potential  is  linear  between  any  two  successive 

points  in  the  potential  table,  the  force  can  be  expressed  simply 

as 

Force  =  (U(X)  -  U(X  +  A)) -(const) 

where    the   constant    converts    the   force   to  proper    units. 

The    determination    of   a    force    or    potential    from    the   tables    is 

just    a    simple   matter    of    interpolation.       For    example,    if   the   atomic 

separation    between    two  atoms    is    1.275a      one   would    just    make   a 

H 

linear    interpolation    between    the    values    at    1.27a      and   1.28a    . 

H  H 

These  tables  are  used  in  both  SUBROUTINE  STEP  which  calculates 

forces  and  SUBROUTINE  IONPE  which  calculates  potential  between 

the  ion  and  the  atoms . 

Two  points  should  be  noted.   First,  the  nearest  neighbor 

separation  is  approximately  4 -436a  .   Therefore,  for  a  separation 

H 

of   4.434a       the    interpolation    is    between    the    value    in    the    table 
H 

at    4- 43a      and   the   value   at    4. 436a     .       Secondly,    the   average    force 
H  H 

between    two  points    is    assumed   to    occur    halfway    between    them.      The 
interpolator    for    the   forces    accounts    for    this    fact. 

One   major    problem   remains    unsolved    in    the   construction    of    the 
force    table    (not    in    the    interpolation    of    forces    in    SUBROUTINE   STEP 
after    the    table  has    been    constructed).       While    the  potentials 
steadily    increase   as    the    interatomic    distances    decrease,    the 
forces    in    the   force    table   do   not    monotonically    increase   with 
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decreasing    atomic    separation.      This    is    true    only   when   potentials 
derived   from    the    Interaction    Potential    Functions    are   used.      No 
difficulty    arises    in    the   region   where   a    Born-Mayer    type    of    ex- 
ponential   fit    for    the   potentials    is    made.       Here   the   forces    behave 
as    expected. 

The   reason    for    this    discrepancy    is    found    in    the   calculation 
of    the    interaction   potentials    that    are   used   by    the   computer    pro- 
gram  as    data.       The    interaction  potential    is    calculated   exactly 

from   0.01    to   0.20a      and   every    0.20a       to    5.0aTi    (i.e.    0.2,    0.4, 

H  H         H 

0.6a   would  be  calculated  exactly,  0.50a   would  not).  Potentials 
H  H 

in  between  these  values  are  calculated  by  linear  interpolation. 
Hence,  the  forces  in  the  table  would  remain  essentially  constant 
in  the  region  where  the  potentials  were  derived  by  linear  inter- 
polation and  would  only  change  abruptly  when  a  point  is  reached 
where  the  potential  was  calculated  exactly.   Thus  while  the 
overall  force  would  increase  as  interatomic  separation  decreased 
there  would  be  regions  where  the  force  stayed  essentially  con- 
stant or  even  decreased  with  decreasing  atomic  separation. 

Two  main  efforts  were  made  to  eliminate  this  discrepancy  both 
of  which  began  by  considering  only  the  potentials  that  were  cal- 
culated exactly  and  discarding  those  that  were  determined  by 
linear  interpolation.   First,  subroutine  ALI  and  ATSE  from  the 
scientific  subroutine  package  of  IBM  were  used  to  fit  known  po- 
tentials with  a  polynomial  and  using  this  polynomial,  points  in 
between  the  known  potentials  were  interpolated.   In  the  second 
method  an  exponential  fit  of  the  known  potentials  was  made,  and 
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an  interpolator  constructed  as  follows.   Assume  the  potential  P 
at  two  points  x   and  x   are  known  and  are  to  be  joined  by  an 
exponential.   We  desire  to  know  the  potential  at  x  . 


+  X 


Xl    Xa    X2 


X2  Xl  . 

AX  =  — — —  where  N  is  an  integer 

N 

(1)  P1    =    Aexp(-Bx1) 

(2)  Pa    s    Aexp(-B(x1+xa))    or    P&    =    P;L    exp(-Bxa) 
Similar ly 

(3)  P2    =    P1exp(-BNAX)    or    ln(P2/P    )    =    -B-NAX) 

(4)  BAX    =    -    ln(P2/P    ) )/N 

Assume   X      =  M«X  where   M   is    a    constant    0   ^M   ^N    then 
a 

P       =    P      exp( -BAX-M) 

Both   the  polynomial    and   exponential    fit    of    known   potential 

yielded  approximately    the   same   results.       In    general,    for    atomic 

separations    of    less    than    1.6a      the    forces    increase  with   decreasing 

H 
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atomic  separation  as  desired.   For  separations  greater  than  1.6a 

H 

the  forces  act  as  expected,  except  that  in  the  immediate  vicinity 

(~.02a  )  of  a  calculated  potential  the  forces  decrease  slightly 
H 

and  then  increase  again.   Since  there  was  little  to  choose  between 
the  two  methods,  the  exponential  interpolator  was  selected.   The 
interpolator  constructs  the  tables  slightly  faster  (~  2  sees)  and 
requires  less  storage  space. 

The  effects  of  this  problem  on  the  force  tables  are  as 

follows : 

+6       .  +4 

a)  Xe      -Si      -no    effect.       The    forces    behave   as    expected    since 

the    Born-Mayer    approximation    is    constructed   beginning   at    1.2a    . 

H 

+4       +4 

b)  Si      -Si         (used  as    bullet-lattice    interaction   potential 

in   preliminary    studies    of    (100)    channeling).      No    effect,    as    a 

Born-Mayer    approximation    is    constructed   beginning   at    1.2a    . 

H 

c)  Si   -Si  '  (used  as  lattice-lattice  interactions  for  all 

directions  of  channeling).   This  force  table  contains  the  error 

discussed  above  for  interatomic  separations  between  1.4  and 

3.0a  .   However,  this  error  has  negligible  effect  on  calculated 
H 

channeled  ranges  since  (a)  there  is  no  error  in  the  force  table 

at  atomic  separations  of  4»4a  ,  the  approximate  nearest  neighbor 

distance  silicon  atoms  and  (b)  the  bullet  passes  the  lattice  atoms 

so  quickly  that  when  two  lattice  atoms  have  an  interatomic  spacing 

of  less  than  3.0a   the  bullet  is  no  longer  under  their  influence. 

H 

While  this  error  in  the  force  table  has  virtually  no  effect 
on  channeling  range  estimates,  before  the  force  table  concept  is 
used  for  other  computer  studies  this  problem  must  be  solved. 


40 


APPENDIX   C    (PROGRAM   PARAMETERS) 

The   following   parameters    have    been    added   to    the   basic    computer 
channeling   model    (Chan   69MODO). 
BPNT    -    Bullet -Lattice  Potential    Energy 
TPNT    -    Lattice-Lattice   Potential    Energy 
BFORCE    -    Bullet -Lattice   Force 
TFORCE    -    Lattice-Lattice   Force 
KA    -    an    integer    constant   which    identifies    the   point    in    the 

potential    table  where   a    Born-Mayer    potential    of   form, 

-  Bx 
Ae         ,    begins 

BPCNST    -    the   value    of   B    in    the    Born-Mayer    Potential    approximation 

for    the    Bullet -Lattice  Potentials 
TPCNST    -    the   value    of   B    in    the    Born-Mayer    Potential    approximation 

for    the   Lattice-Lattice  Potential 
ICHECK,    ICOUNT,    NCOUNT    -    constants    used    in    the    exponential    inter- 
polator   between    calculated    interactions    potentials    in 
constructing    the   Potential   Tables 
BDELX   -    a    variable    (BAX)    (see   appendix    II)    used    in    the   exponential 
interpolator    between    interaction   potentials    in    con- 
structing   the   Potential    Table 
ITITLE    -    80   alpha    numeric    characters    used    on    a    heading    card   for 

potential    data    deck 
LCUT    -    integer    constant   with    value   of   either    0    or    1 .       Used   to 
determine   which   atoms    are    in    the   core    (0)    or    shell    (1) 
of   the   crystal 
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CVEDA    -    conversion    factor    used   to    convert    electron    volts   per 

Bohr    to   newtons   per    meter 

CLU    -    conversion    factor    used    to    convert    distances    in    LU    to 

distances    in    units    of    100a 

H 

AFC    -    distance    in    100a 

H 

IFC  -  integer  value  of  AFC 

CLU,  AFC,  and  IFC  are  used  in  the  following  way.   Assume 

interatomic  spacing  between  atoms  =  0 . 5LU 

AFC  =  (CLU)  x  (DIST)  =  2  56 .14 

IFC  =  256 

An   AFC    of   256.14    corresponds    to   2.56l4a 

H 

An    IFC    of    256    corresponds    to   2.56a 

H 

AFC1    -    nearest    neighbor    distance    in    units    of   100a 

H 

AFC2    -    the  reciprocal    of    the    difference   between    the   nearest 

neighbor    distance   and   443. 0 «       Used    in    Subroutines    STEP 

and    IONPE. 

DFF1  -  the  difference  in  100a   between  the  atomic  separations  of 

H 

the  bullet  and  lattice  atom  and  the  nearest  neighbor 
distance  of  silicon 
GETIME,  SETIME  -  machine  subroutines  used  to  determine  actual 

program  running  time.   They  terminate  program  execution  and 
have  cards  punched  before  exceeding  the  running  time  allo- 
cated to  the  program  by  the  computer.   This  allows  a  re- 
start capability. 
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Computer    Program    -    SUBROUTINE    DL100 

*******************       SUBRPtlTINF    0L10O       ********************* 

*  THIS     IS    A     LATTICE    GENERATOR    FOR    THf=    DIAMONO    LATTICF     IN       * 

*  THF     (IOC)     rRIcNTATIPN.  * 

*  THE    CRYSTAL    PcVcl_PPS    IN    T"c    OPOER    X,FOLLOWEP    BY     7,  * 

*  FPLLOWFR    PY    Y  * 
************************************************************ 

COMMON /C 0*1 /RX( 200)  ,"Y(2  00)  ,R?(200  ) ,LCUT(200  )  .1  L, LC 

fPMMrN/CCV2/ROF  ,RCF2  ,ROFV,  AC,PAC,PPT,PTCt  RFPTr  ,  F^TC, 
2    FM,pPIV?TO0T 

CHMMPN/f  TM4/IX,  I  YtI7flXP,lYPtI7D,sCX,SCYfSC7tIl>rcO 

COMMCN/CCM7/Rl,LSS,SPX,S°7tCOX,COY,C07 

COMMpN/CC^/     ICX,IOYtIO7(2  00),nxT,nYTtr7TtTDCX,Tf,0Yt 
?    TP07 

COMMON  /cr*?_2  /XLL,Xl«r,7LL,7LF 
************************************************************ 

*  tmp  FHLLTWINC  rGNSTANTS  ARE  USEP  FITHF9  BY  THIS  SMP-     * 

*  pntjTTNF    TO    GENERATE    THF    LATTICE    OR    IN    OTHFP    SUBROUTINES    * 

*  4  S    "FQIIIPEP.  * 

************************************************************ 

XLL  =  C.<=o 
XLF=3.C1 
7LL  =  CS<= 
7\  F=3.oi 
M  =  2 
L6=2 
IB=l 

iox=i 
ipy=m 

TX  =  <5 

IY=5 

I7  =  * 

CCX=1.C 

SC  Y=  1  .  C 

SC7=1.0 

nXT=?.C*cCX 

DYT=2.C*CCY 

PZT=2.0*SC7 

0Y=0.5 

PX=0.5 

07=0.5 

*********************************************$********•<***** 

*  THF    cPLLOWING    SERIFS    OF    NESTEO    PC    LOOPS    GENFRATE    THC  * 

*  LATTICE    ANP     IN^LRF     IT    CONTAINS    THE    Rc0UIRE0    SY^CT"Y  * 

************************************************************ 

JT  =  0 


Y=-SCY 

PP    6C    J=1,IY 

Y=Y+SC  Y 

KT=0 

7=-SCZ 

DO    50    K=1,I7 

Z=Z+SC7 

1^=0 

X=-SCX 

DP    AC    !=l,TX 

x=x+scx 

ITT=IT>JT>KT 

IF(ITT-CITT/2*2)1    40,30,40 

30 

PXI  M)  =  X 

o Y( M)=Y 

R7( M)=7 

ID7(M|=5 

L  =  M 

M  =  M  +  1 

40 

IT=IT+l 

50 

KT  =  KT4-1 

IF(   I.E0.5)    PC    tp    80 

On    200    IA=LA,L 
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202 
203 


200 


PC 


200,200,202 
200,200,203 


TF(  4.C*SCX-RX(TA)*SCX) 
TP(  ^.C*5T7-R7(  IA)  *SCZ) 
o  X(l+IP)=RX(  iA)*nx 
PYI L  +  TP )=9Y( IA )+DY 
R7(  L+IR)=P7(  IAJ+07 
1^7 ( l+TP)=4 
IB=IP+1 
C^NTINI!^ 
u  =  L  +  <5 

LA=^ 
TR=1 
JT=JT+I 
tl=M-l 

*  THE    pOlinwTNG    Of    I  0D&    DIVinFS    -he    CRYSTAL     I  NT  H    AN    ACTIVE* 

*  CORP    AND    AM     INACTIVE     SURROUNDING    ^H«=LL.  * 

A*************^*^^**^******^*^*** *************************** 

nn   2   i=2,ll 


OR.RX(T) .L 

LCtJTfT  )=1 

nR.R7(I  )  .ft 


LCIIT(  I  )=C 

Tp(PX(  I  LOT.  3.1. 

TM&Yf  T  )  .HT.2.5) 
IC(P7(  I  ).li.C9. 
CnNTIkWE 
LCUT< 1)=C 

Rl  =  [ppr+0Y(2}|  /APS (COY) 
t?ETl!RN 

CMp 


,Q   ) 
I    ) 


LCUT(  I  )  =  l 
LOJT( T)=l 
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Computer    Program    -    SUBROUTINE    DL110 

*******************        SURRPUTINF     nLllO        ********************* 

*  THIS    !S     ft     LATTICE    OENFRATOR    FOR    THE    DIAMOND    LATTICF     IN       * 

*  THF     (110)     CRICNTATI ON.  * 

*  THF  CRY$Tflt  0<=\/,-LnPS  IN  THE  CR0CR  X,C0L10WFP  RY  7,  * 
*PnLLnWFHP,  YY  * 
*************************************************** ********* 

fOMMPN/rrvi/Rxi 2  00)  ,RY(200)  ♦  P 7 { 200 ) , LCUT  I  2 30  )  ,  LL , L 0 
CnMMpNj/rrv?/PrP  TRpP2  ,RQEVf  AC  ,PAC  ?PFTC,  PTft  DFPT C  ,  FPTC, 
2    FM.ocj v.T-OT 
GOMMnN/CpMA/IX,  I  Y,I7tIXP,IYP,I7-2fsCX,SCY,SC7,IPFFP 

cn^MPM/crM7/R!  fLSS,SPx,sDz,cox?coY,cp7 

COMMPN/f  CV<V    IPX,I0Y,IP7 (200)  ,  0  XT  t  PYT  f  0  7^  f  T  or*  ,  j  pnYt 
2    TPp7 
CPMMPN/cry??/Xt.L  ,XLC  ,7LL  ,71  «= 
************************************************************ 

*  THC     FOLIOWIKT,    rr-N$T4NT<     ARE    U^FO     CITHE»     RY    THI*     S'JR-  * 

*  PPUTINF     TP    GENE  PA  TF     THF     LATTICE     OR     IN    CTHFR    9 '  IRPOOT  tnjf«-     * 

SfSPFOHTOFD.  * 

************************************************************ 

XLL=C.*C7 
XLF  =  2.  131 

7LI  =C.c 
7LF=3.Cl 

M  =  2 

Tnx=l 

IOY=23 

IX- 5 

IY=«5 

I7  =  c 

Pa=1.0/SQ«TC2.0) 

SCX=RO 

<C Y=RP 

SC7*1.0 

PXT  =  2. c*cr  X 

0YT=2.C*?C Y 

r>7T=2.C*^C  7 

OX=PP 

0Y=0.C 

07=0. c 
************************************************************ 

*  THE    FOLLOWING    SERIFS    pe    NFSTFn    no    LOOPS    GENERATE    THE  * 

*  LATTTfp     4  ND    INSIJRF     IT    CPNTAINS    thp    PEOUIRFD    SYMMETRY  * 

*********************************** **************** ********* 
)T=0 

Y=- SC  Y 

PO    40    J=1,TY 

Y=Y+SC  Y 

KT=C 

7  =-  SC  7 

on  «sc  k  =  i,i7 

7=7+SC7 

IT  =  0 

x=-sr  x 

MA=M 

00  ac    1  =  1, IX 

X=  X+  CC  X 

I*( Tf-( IT/2)*2)  21,11,21 

11  TF( JT-( JT/2) *2)  40,12,40 

12  IF(KT-(KT/2)*2)  ^o, 30,40 

21  IF( JT-( JT/2) *2)     22,^0,22 

22  IF(KT-(KT/2) *2)     30,40,30 
30    PX(  ^)  =  X 

P Y( M)=Y 
R7( M)=7 
L  =  m 
V  =  M+l 
AO     IT=IT-H 

MR=M-1 
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202 
2C3 

POO 
21* 


710 

*o 


2 

***** 

*  TH 

*  fr 
***** 


Ic( 
Ic( 
!F( 

nn 

IP( 

TF( 
"X( 
RY( 
P7< 

M  =  M 

rr»N 
r,n 

LP  = 
PX( 

5Y( 
P7( 

M  =  M 

GO 
COM 
KT  = 
JT= 
LL  = 

HP, 

107 
**** 

F  P^ 
QF  /» 
**** 

no 
Lru 

IF( 
IF( 
IF( 
CON 
LCD 
Rl= 

CMP 


MR- 
JT. 
JT. 
20C 
3.c 
^.0 
M  )  = 

V)  = 
M)  = 
+  1 

TIN 

TO 

M-2 

M)  = 
M)  = 

M)  = 

+  1 

Tp 
TIN' 
KT4- 
JT  + 
M-l 
2    I 

( I ) 

*** 
LL° 

Njp 
*** 

o  T 
T(  I 
p  Y( 
D7( 
»7( 
TIN 

T(  1 

(  co 

II'  \' 


ma  )    21C#210f198 

c0.1)    00    TH    2lc 
F0.3)    GO    TO 
I=MA,MB 

*sr  x-ox(! ) ) 
*cc  7-pJ<  i ) ) 
px(  i  n-nx 

FY( I )+0Y 
P7(  I  )+P7 

UF 

210 

CXI  LP) -OX 
PYC LR) +DY 
P7( L')+D7 


I  cc 
UF 
1 
I 

=  2,1 

_  c 

**** 

WT  NO 
AN'  I 
**** 

=  2,L 
)=C 
I  ).G 

I  ).G 
!  ).L 


215 

200,200,202 
200,200,203 


******$*******$******  ************** ********* 

^    l°HP    PIVIDFS    THE    rPYSTAL     into    &N    ACTIVE* 

N*r.TTVP     ^'JPPQUN^ING    rHcU.  * 

******************************************** 

L 

T.2.2)  LCHT(I)=1 
T.3.^)  LCIJT(I)=1 
t.0.4)     LCUT(I)=1 


IT 

)=C 

C+PY(2) ) /£Bc(rrY) 
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Computer    Program    -    SUBROUTINE    DLlll 

*******************       SURROUTINF    DLlll        ********************* 

*  THIS     TS    A     IATTICE    GENERATOR    FDR    THF    DIAMpNO    LATTICF     IN       * 

*  THF     (  111 )     CRICNTATICN.  * 

*  THF    CRYSTAL    DEVELOPS    IN    THE    ORDER    X, FOLLOWED    BY     Z,  * 

*  FOLLOWEP  RY  Y  * 
************************************************************ 

CPMMniM/CCMl/QX(20O)  ,PY(2  00)  .R7(200)  ,LCUT  (200),LL,LD 
C0MMCN/CCM2/R0E,R0E2 ,  RPF Mf AC  ,  PAT , PPTC , PT C , PFPTC, FPTC, 

?    PMfcpy v.tppt 
COMMpv /ffM&/ 1 X,IY,I7,IXP,IYP,I7P,SCX,SCY,SC7,IDEEP 
COMMnM/rCM7/RlfLSS,SPXfSRZ,COX,COY,COZ 
COMMpM/rrwq/    lDXtInY, 107(200) ,QXT , DYT , D7T , T PCX, T POY, 

2    TPP7 
COMMON /CCM22/ XL L ,XLF t 7 LL  ,7LF , YL L ,Y LF 
************************************************************ 

*  TH«=    FOLLOWING    CONSTANTS    ARF    U^ED    PITHER    BY    THIS     SUB-  * 

*  POUTINE  TC  OENERATF  THE  LATTICE  OR  IN  0THER  S'fBRQW  IMPS  * 
♦ASRFOUIRFO.  * 
************************************************************ 

XLL=0.6°7 

XLF=2.131 

YLL=2.5 

YLF=^.4 

7LL=1.21E 

ZLF=3.^F4 

m  =  2 

L*  =  2 

I«  =  l 

mx=i 

IDY=** 

00    2    1=2,150 
2    ID7( T )=* 
IX=5 
TY=7 

IZ  =  * 

SCX=1.C/SQRT(2.0) 
SCY=2.C/S0RT(3.0) 
SC7=SQPT( 1.5J 
SSC7  =  SC7/3.0 

nxT=2.r*c;c  x 

DYT=2.C*SORT(3.0) 
07.^=2. C*<:CZ 
DX=0.C 
r>Y=SORT(  C.75) 

07=0.0 
*************************************************** ********* 

*  THE    POLLOWTNG    SERIES    0^    NESTEO    00    LOOPS    GENERATE    THE  * 

*  L&TTTCE  A\o  INSURE  IT  CONTAINS  THC  REQUIRED  SY"VETRY  * 
*************************************************** ********* 

JT=0 

y=-  sr  Y 

DP    *0    J=1,IY 

Y=Y+SC  V 

KT  =  0 

jTS  =  JT4-JT/3 

Z=-^C7 

00     ^0     K=1,I7 
IF(K.EO.c)    OC    to    80 
°0    COMTINUF 
Z=7+SC7 
IT  =  0 
X=- SC  X 
DO    40     I=lfTX 

x=x+sc  x 

IN=IT+JT*+KT 

IF( TN-f I*/2)*2)     40,30,40 
30    RX(M)  =  X 
RY(M)=Y 
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33 


PC 


2  0 
2C 


20C 


c 
IC 

*  T 

*  C 


JT 
JT 

IP 
JT 

7P 
R7 

go 

R7 
L  = 

M= 

00 

IF 
IF 
IP 

00 
IT 
KT 
IF 

no 
IP 

IF 

RX 
RY 
P7 
IB 

MA 

CO 

M  = 

LA 
IB 
JT 
LL 


JT-3*(JT/3))     31,35,31 

=  JT 

=JTT-3 

JTT)     33,35,32 

=  JTT4-3 

JTT 

M)=7+7P*$SCZ 
Tn    37 
M)=7 


HF    F 

OP  c 

**** 
nn 
LC 
IP 
IP 
1^ 

3  rn 
LC 
PI 
Pc 
FN 


+  1 

TO 
J. 
J. 
J. 

Tn 

IT 

KT 

J. 

0 

20 

4. 

A. 

L  + 

L  + 

L  + 

IP 

M.A 

TI 

+  M 

M 
1 
JT 

LL 

NO 

3 

T( 

P  X 
P7 
PY 
TI 
T( 
(P 
IJP 


AC 
c0.1) 

FQ.4) 
F0.7) 

P5 
+  1 

+  1 

c0.7) 


OG 
GC 
GC 


TO 

Tn 

TD 


90 


GC    T0    100 


0     T  J  -14     I 

c*cr  x-p  x(  I.*  )  *SCX)    200,230,202 

C*?C7-D7( i a ) ) 200,203,203 

IP)=RX( IA)+0X 

IP )=PY( IA )+DY 

I8)=P7( IA)+07 

♦  1 

+  1 

NLF 

\+l 


+  1 
I 

*** 
PVT 
AM 
*** 
1=2 
I  )  = 
(  T  ) 
{  T  ) 
(  T  ) 
NLF 
1)  = 
Oc  + 
N 


NG    OC    I OCP    DTVIHFS    THE    CRYSTAL    INTO    AN    ACTIVE* 
INACTIVE    SURROUNDING    ^HFLL.  * 

********  ***************  *****#**##*#*  ********** 

tLL 


.LT.C.7.CR.'X(I ) .GT. 2.2) 
.LT.1.4.PR.R7(I ) .GT.3.^) 
.GT.5.2)     LCUT(I)=1 


PY<2) ) /ABSCCHY) 


LCUT(I)=1 
LCU^( I )=1 
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Computer    Program    -    Exponential    Interpolator 


****** 

*  PEA 

*  MAY 
****** 


4C5 

****** 

*  CON 

****** 

410 

415 

****** 

*  MJK 

*  pnr 
****** 


416 
417 


41$ 

421 

****** 

*  PET 

*  POT 

****** 

422 


425 

420 

****** 

*  CON 

****** 

430 

****** 

*  PCA 

*  MAY 
****** 


445 

****** 

*  CON 
****** 

450 

4C  5 
****** 

*  MAK 

*  POT 
****** 

456 
457 


***** 
0  LAT 
EP  ftp 
***** 
RFAP( 
PF*0( 
RFAO( 
PFA0( 
IF  (K.I 
ICHFf 
r,o  t 
***** 

STPur 

***** 
DP  41 
TPNT( 
***** 

F    6N 

***** 

IcmiN 
no   m 

ncpun 

BOFLX 
TDNJT( 

IF(  IC 

icnuN 

OP  to 
***** 
FRMIM 
FNTU 
***** 

TTPNT 
00  ^2 
TPNTf 

IF(  TO 
TPMT( 

CONTI 
***** 

STRUC 
***** 

DO    43 

TFORC 
***** 

D  RUL 
FR  AP 
***** 
REAP( 
RFftO( 
PPAD( 
RFAO( 
IC(KA 
ICHEC 
GO  TO 
***** 

STRUC 
***** 
DO  45 
BPMT( 
***** 
F  ftN 
FNTIA 
***** 
ICHFC 
10  PUN 
00  45 
NC  0  UN 


***** 

TlfF-LA 
DRfXI Vfi 
******* 

5,9°2S) 
5,c«10) 
^,CP15) 
5,CP2C) 
)  405,4 
K=4*C 
416 

******* 

T  A  ROR 
******* 

c  I=K4, 
1  )=TPM 
******* 

CXPCN=N 

LS 

******* 

K=KA 
T  =  2C 
c    1=1,2 

T  =  ir cun 

=-ALCG( 
TCPUNT+ 
HFfK-NC 
T  =  NCPLN 

417 
******* 

F     PPTCN 

L  TABLE 
******* 

=TFNT(  I 
C  1=1,5 
I )  =  ( TON 
NT< I  )  )4 

t  )=c: 

NUE 
******* 

T  A  LAT 
******* 
0     1=1,4 

c( T)=(T 
******* 

LFT-LAT 

DPPXIMA 

******* 
5,<=825) 
5,9*10) 
5,cei*) 
5,CR2C) 
)  445,4 
K=46C 

456 
******* 

t  A  BOP 
******* 

5  T=KA, 
! )=BPNT 
******* 

EXPONEM 

LS 
******* 

K=KA 

T  =  2C 

<3     1  =  1,2 

T=ICCLN 


FXPnNFNTIAL  INTERPOLATOR  **************** 
TTTfF  POTENTIALS  AND  DETERMINE  Ic  A  BORN-* 
TinN    I*    NECESSARY  * 

****************************************** 

(ITI TLF(I)  ,1=1 ,20) 
(TPNT(I)  ,1=1  ,20) 
(TPNT( I) ,1=21 ,500) 
KA,TPCNST 
C5,410 


****************************************** 

N-MAYER    DCTENTIAL  * 

****************************************** 

500 

(KA)  *FXP(TPCNST*(  ( I/100.)-(KA/10O.)  )  ) 

****************************************** 

TIAL    INTERPOLATION    R«=TWFFN    CALCULATFC  * 

* 
****************************************** 


:c 

IT+20 

TPMT(NCOONT) /TPNTf I  COUNT ) )/2C . 

I )=TPNT( ICOUNT)*EXD( (-I)*BDELX ) 

GliNT)    422,422,421 
IT 


****************************************** 
TIAL  AT  NEAREST  NEIGHBOR  DIST ANCF , 7FP0  * 
,AN0  CONVERT  POTENTIAL  TO  ELFCTRnN  VOLTS  * 
****************************************** 

CC)-M(TPNT(IFC  +  1)-T?NT(  IFC)  )*(  AFC-IFC)  ) 

00 

T(  I  )-TTPNT)*i3.^ 

25,420,420 


****************************************** 

tTCE-LATtice    FORCE    TABLF  * 

****************************************** 

99 

PNT(  I  )-TP\T(H-l)  )*CVFOA 

****************************************** 

TICE     PCTFNTIALS     AND    PETERMINF    IF    A    BORN-    * 
TION    I  <=    NECC<:SARY  * 

****************************************** 

(ITITLE(I)  ,1=21,40) 
(BPNT(I) ,1=1 ,20) 

(BPNT(I) ,1=21,500) 

KA,BOCNST 

45,450 


****************************************** 

N-MAYER    POTENTIAL  * 

****************************************** 

500 

IKA)*FX°<BDCNST*( (1/100. )-(KA/100.) ) ) 

****************************************** 

TIAL    INTERPOLATION    PETWEEN    CALCULATED  * 

* 
****************************************** 


0 
T+2  0 
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4f  1 


**** 

*  D 

*  P 
**** 

46 


^5Q     RP 

TF 
TF 
p.n 

**** 

FTCD 

OT=M 

**** 


47 

4P 

**** 

*       C 

**** 

48 

4C 
**** 


2    RR 

nn 

RP 

IF 

C      RP 

C    CP 
**** 

**** 

c  bf 
**** 


FLX=-H 
NT(  ICC1! 

(  TCHFCK 

ni)NT=\r 

Tp    i  =  7 
******* 

M  T  M  c     F  P 

******* 
PNT^PFNi 
^80  !  = 
NT| I )  =  { 
IRPMTII 
Nt( I )=C 
NT!  NIJF 
******* 

r  i  ir  T  a 

******* 

4CQ     1  = 

P"CE(  I) 
******* 


n^f RPNT 

-NCCUNT 
HUNT 

******* 

TcNTf  A l 

RLC  ,ANP 
******* 

TC  IFCK 

1,500 

*0NT( I) 
>)475,4 


(NCOUNT)  /RPNTU 
°MT( ICOUNT) *EXC 
)    462,462,461 


*************** 
AT  NE^RF^T  NEI 
PPNVFRT    PATENT 

*************** 

( (B°NT(IFOll-B 

-BBPMT)*13.6 
BO, 480 


******* 

BULLET- 
******* 

4PQ 

RCNT( 
END 


U 

** 


*************** 

LATTICE    ^CPCE    T 
*************** 

I)-BPNT(H-1  )  )*F 

PIF     INTF.P.POI  ftTH 


CDJMT))/20. 
f (-I)*BOFLX) 


*********************** 
PHBOR    DISTANCE, 7ERO  * 

ML    TQ    ELFCTP0N    VOLTS    * 
*********************** 

PNT( IFC) )*(AFC-IFC) ) 


*********************** 

*BLE  * 

*********************** 

VEDA 

o  ******************** 
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Figure   2 
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Figure  3 


53 


<lll> 


<l  00> 


Figure  4.   The  shape  and  sizes  of  the  major 
channels  in  silicon. 
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